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Abstract  We  proposed  a  technique  for  reconstructing 
bio-electric  conductivity  distribution  from  measured  tangential 
magnet  field  data.  That  is,  in  the  2D  ROI  with  an  arbitrary 
depth  two  orthogonal  tangential  components  of  3D  current 
vector  field  were  respectively  estimated  from  the  synthetically 
measured  two  tangential  magnetic  fields,  from  which  the  2D 
conductivity  distribution  was  estimated.  To  cope  with  inevitable 
measurement  errors,  a  robust  numerical  solution  was 
developed,  in  which  the  mollification  method  and  the 
regularization  method  were  efficiently  utilized.  By  performing 
this  2D  reconstruction  at  each  depth,  3D  conductivity 
distribution  could  be  estimated.  As  typical  applications, 
pathological  state  and/or  bio-electric  conductive  path  can  be 
evaluated  by  monitoring  the  temporal  change  of  the 
reconstructed  3D  conductivity  distribution.  The  feasibility  of 
this  technique  was  briefly  verified  by  reconstructing  a 
conductivity  distribution  using  simulated  noise-filled  magnetic 
field  data,  with  resultant  reconstructions  indicating  that  the 
approach  could  provide  a  practical  means  for  robustly 
reconstructing  a  conductivity  distribution. 
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I.  Introduction 

Utilizing  SQUID  magnetometer,  ferromagnetic  particles  or 
electrical  activities  are  evaluated  (e.g.,  dust  in  lung  [1], 
magnetocardiogram  (MCG)  [1,  2],  magnetomyogram  (MMG)  [3], 
and  magnetenchepharogram  (MEG)  [1]:  visually  evoked  [4], 
auditory  evoked  [5],  somatosensory  evoked  [6],  spontaneous 
motion-related  [7],  cognition-related  [8],  and  electromagnetically 
excited  neuron-  and  nerve-induced  [9]).  Often,  by  modeling  tissue 
as  a  conductor,  current  sources  [10,  11,  12,  13]  or  current  fields  [2,  9, 
14,  15,  16-18]  are  estimated.  Various  current  field  estimation 
techniques  are  developed  (e.g.,  the  unique  minimum-norm 
least-squares  estimation  method  [14,  18].  the  Wiener  estimation 
method  [15],  the  Fourier’s  filtering-based  method  [16,  17]).  From 
the  measured  normal  magnetic  field  to  the  reconstruction  plane,  any 


current  field  components  cannot  be  determined  [2,  9,  14,  16].  Then, 
the  tangential  magneticometer  [3,  19,  20,  21]  is  utilized  [17,  18], 
That  is,  utilization  of  the  magnetic  fields  of  orthogonal  two 
tangential  components  realizes  unique  determination  of  the  2D 
current  vector  field. 

Recently,  on  3D  object  the  synthetic  aperture  technique 
is  applied  to  the  normal  magnetic  field  measurement  [22],  Thus, 
novel  application  of  the  synthetic  aperture  technique  to  the 
tangential  magnetic  field  measurements  could  allow  determination 
of  the  3D  fields  of  orthogonal  two  tangential  components  of  3D 
current  vector  field.  In  this  report,  to  allow  providing  the 
quantitative  information  about  the  pathological  state  and/or 
bio-electric  conductive  path  (e.g.,  neural  networks  etc),  we  present  a 
reconstruction  technique  for  monitoring  the  temporal  change  of  the 
3D  conductivity  distribution  using  only  the  estimated  tangential 
current  field  data.  We  previously  proposed  a  differential  type 
conductivity  inverse  problem  using  either  current  or  potential  field 
data  [23-25],  in  which,  if  no  current  sources  exist  in  the  region  of 
interest  (ROI),  with  respect  to  reference  conductivity  values 
provided  in  an  arbitrary  point  [23,  24]  or  regions  [25]  the  spatial 
distribution  globally  of  the  relative  values  can  be  determined  from 
measured  sufficient  number  of  independent  fields  of  another 
physical  variable  in  that  region.  Thus,  compared  with  an  eddy 
current  model  [26]  previously  proposed  in  the  field  of  NDE,  in 
principle  our  technique  can  deal  with  an  arbitrary  geometry  of  the 
target  and  arbitrary  current  sources  existing  outside  ROI. 
Furthermore,  we  delineated  that  proper  configurations  of  the 
reference  regions  ami  the  externally  situated  current  sources  allow 
reducing  the  required  number  of  independent  current  field 
measurements  for  uniquely  determining  the  conductivity 
distribution  [25],  This  approach  is  particularly  effective  in  this  case 
when  an  inherent  current  flow  exists  there  and  a  fewer  field 
measurements  are  required.  In  real-world  applications,  however,  the 
problem  inevitably  becomes  ill-posed  due  to  inevitable 
measurement  errors  [24]  and  improper  configurations  [45]  making 
it  impossible  to  guarantee  the  existence  of  a  stable  and  unique  target 
global  distribution.  However,  the  demonstrated  ability  using 
simulated  noise-filled  current  field  data  under  the  improper 
configurations  of  our  developed  numerical-based 
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implicit-integration  approach  to  reconstruct  the  globally  relative 
conductivity  distribution  is  concluded  to  confirm  the  potential  value 
that  it  can  yield  the  acceptable  unique  reconstruction  even  when 
fewer  field  measurements  are  unfortunately  taken  [44],  To  cope 
with  inevitable  measurement  errors,  in  conjunction  with  the 
synthetic  aperture  technique  a  robust  current  field  estimation 
method  was  developed,  in  which  the  mollification  method  and  the 
regularization  method  were  efficiently  utilized.  By  enhancing  and 
combining  the  current  estimation  method  and  the  conductivity 
estimation  method,  we  obtained  the  3-D  conductivity  reconstruction 
technique.  The  feasibility  of  the  technique  is  evaluated  by 
reconstructing  a  globally  relative  conductivity  distribution  using 
simulated  noise-filled  magnetic  data,  with  resultant  reconstructions 
indicating  that  the  approach  provides  a  practical  means  for  robustly 
reconstructing  conductivity  distribution. 

II.  Conductivity  reconstruction  technique 

The  conductive  target  medium  existing  in  3D  space  expressed  by 
Cartesian  coordinate  system  ( x ,  y,  z)  has  3D  conductivity 
distribution  G  (x,  y,  z).  The  3D  ROI  £2  is  taken  into  the  target 
medium.  The  conductivity  distribution  G  (,v,  y,  z)  in  £2  is 
spatially  smooth,  and  the  reference  conductivity  values  G  ,  (1=1, 
...  ,  N  )  aie  provided  in  reference  regions  CO  ,  existing  in  £2  . 
i.e., 

In  G  (x,  y,  z)  =  In  G  , , 

(x,y,z)  <=  CO,  (/=  1 . JV).  (1) 

Under  the  conditions,  magnetic  field  B  (x,  y,  z)  (=  (Bx  ( x ,  y,  z),  Byfr, 
y,  z),  Bz  (,v,  y,  z))T)  is  generated  by  3D  current  density  field  J  ( x,y,z ) 
(=  (Jx  (x,  y,  z),  Jy  ( x,y,z ),  Jz  (x,  y,  z))T)  due  to  arbitrary  current 
sources  existing  inside  and/or  outside  the  ROI  £2  .  Here,  we 
assume  that  the  following  first-order  PDEs  hold  in  £2  for  the  3D 
quasi-relative  conductivity  distribution  In  G  (x,  y,  z)  [23],  i.e., 

V x  J(x,y,z) 

+  Vln(l/  g(x,  y,z))x  J(x,  y,z)  =0.  (2) 

If  synthetic  measurements  could  be  taken  at  the  infinite 
plane  of  [z  =  Zo  ,  —  00  5s  jc  5s  oo>— 00  <y<  00} 
over  the  medium  without  no  errors  on  the  magnetic  field 
components  Bx  (,v,  y,  zy,  Z)  and  By  (x,  y,  zy,  Z)  respectively 
generated  by  the  current  density  field  components  Jy  (x,  y,  Z)  and  Jx 
(x,  y,  Z)  in  an  arbitrary  plane  of  z  =  Z  in  the  target  medium,  by 
mathematically  resolving  the  Biot-Savart's  law  [16,  17],  i.e., 

B  (x,y,zy,Z) 

-j-  JU\/,z')x— clS  ,(3) 

An  \  r  I3 


\ 

Z'=Z,-o°<x‘<oo,-oo<y‘<oo 


the  2D  fields  of  the  components  Jx  (x,  y,  Z)  and  Jy  (x,  y,  Z)  of  the 
3D  flow  can  be  uniquely  determined,  where  permeability  /O  is 
the  given  spatial  function,  and  r  is  the  distance  vector  from  (x\  y\ 
Z)  to  (x,  y,  zo)-  That  is,  in  this  ideal  case,  by  performing  the  Fourier's 
analysis  on  every  plane  the  3D  fields  of  the  components  Jx  (.r,  y,  z) 
and  Jy  (.r,  y,  z)  can  be  uniquely  determined.  Thus,  in  practice, 
coordinate  system  should  be  taken  such  that  generated  current  flow 
is  dominant  in  the  estimation  plane. 

To  determine  the  3D  conductivity  distribution  G  (x,  y, 
z),  one  of  three  PDEs  (2)  is  handled  for  every  plane,  i.e.,  for  the 
plane  of  z  =  Z, 


(jy  (x,y,Z) 


-Jx  (x,y,  Z)) 


(d 


ln<r(x,  y,Z) 
lncr(x,  y,Z) 


J 


=  ^-Jy(x,y,Z)  +  ^-Jx(x,y,Z).  (4) 
ox  ay 

When  in  the  plane  only  single  current  field  of  components  Jx  (x,  y, 
Z)  and  Jy  (r,  y,  Z)  can  be  measured,  to  allow  unique  determination 
of  the  2D  conductivity  distribution  In  G  (x,  y,  Z),  the  reference 
regions  crossing  (0 ,  (l  =  1,  ...  ,  N)  must  extend  in  the  parallel 
direction  of  the  imitated  2D  current  flow  (Jx  (x,  y,  Z),  Jy  (x,  y,  Z)) 1 
[25].  Fortunately,  if  measurements  of  two  independent  current  fields 
Ji(x,  y,  z)  and  J2  (.r,  y,  z)  can  be  taken,  correspondingly  derived  two 
PDEs  are  handled  as  simultaneous  equations  [23,  24],  In  this  case, 
existence  in  the  plane  of  only  a  reference  point  enables  the  unique 
determination  [23]. 

In  practice,  however,  to  determine  the  current  field 
components  Jx  (x,  y,  Z)  and  Jy  (x,  y,  Z)  in  the  finite  2D  ROI  £2  z, 
i.e.,  the  finite  2D  region  of  z  =  Z  in  the  3D  ROI  £2  .  measurements 
of  the  magnetic  field  components  Bx  (x,  y,  zy,  Z)  and  By  (x,  y,  zy,  Z) 
are  performed  in  the  finite  2D  region  of  z  =  Zo  over  the  discrete  2D 
ROI  £2  K  on  a  discrete  Cartesian  coordinate  system  (i,  j,  K)  in 
which  i,  j ,  and  K  are  arbitrary  integers  respectively  satisfying  (,v,  y, 
Z)  =  (i  A  x,  j  Ay,  K  A  z)  with  respect  to  the  sufficiently  small 
spatial  increments  A  x,  A  y,  and  A  z  in  the  x-,  y-,  and 
z-directions.  As  the  measured  magnetic  field  data  Bx  (i,j,  ky,  K)  and 
By  (i,j,  ky,  K)  (k0  =  Zo/  A z)  are  contaminated  with  measurement 
errors,  the  current  field  components  Jx  (i,  j,  K)  and  Jy  (;,  j,  K)  are 
estimated  in  the  numerical  least  squares  sense.  Here,  one  should 
note  that  the  current  sources/current  flow  outside  the  ROI  £2 
significantly  contributes  to  magnetic  fields  over  the  neighborhood 
of  the  border  of  £2  ;  then  resulting  in  the  estimation  errors  in  J  (i,j, 
K).  Thus,  the  extra  contributions  in  Bx  (i,j,  k0‘,  K)  and  By  (i,j,  ky,  K) 
are  reduced  by  being  windowed  with  the  weight  function  W (i,j,  ky. 


K ),  i.e.,  the  spatial  function  having  absolute  values  of  the  magnetic 
field  B  (/,  j,  kg,  K)  calculated  from  uniform  2D  current  flow  with 
unity  magnitude  in  an  arbitrary  direction  in  the  2D  ROI  £2  K.  In 
addition  to  the  fact  that  the  current  field  components  Jx  ( / ,  j,  K)  and 
Jy  0,  j ,  K)  and  their  2D  gradients  are  bounded,  utilization  of  the 
regularization  method  under  the  assumption  of  the  boundness  of 
their  2D  Laplacians  enables  stable  estimation  of  the  3D  fields  of  Jx 
O’,  / ,  k)  and  Jy  ( i ,  j,  k).  That  is,  minimization  is  performed  in  every 
2D  ROI  £2  K  of  the  following  functional  with  respect  to  the  raw 
vector  J  comprised  of  both  the  unknown  current  densities  Jx  (i,j,  K) 
and  Jy  (i,j,  K),  i.e., 

Rj  (j)  =  iibw  -  ljii2  +  a  j,  iijii2  +  a  J2  iidjii2  +  a  „  iid'djii2,  (5) 

where  the  raw  vector  Bw  is  comprised  of  the  lexicographically 
ordered  windowed  magnetic  data  W (i,j,  k0 )  Bx  ( i,j ,  k0;  K)  and  WO, 
j,  k0)  By  0,  J,  kg,  K),  the  matrix  L  is  of  the  lead  field  data,  the  matrix 
D  is  of  coefficients  of  employed  first-order  approximations  of  2D 
partial  derivatives  of  Jx  ( / ,  j,  K)  and  Jy  0,  j,  K ),  and  OC  ji,  OC  J2,  and 
OC  J3  are  so-called  regularization  parameters  adjusting  the  relative 
weights  for  Rj  (J)  of  the  correspondingly  derived  penalty  terms,  i.e., 
L2  norms  of  J,  of  the  2D  gradient,  and  of  the  2D  Laplacian. 
Providing  that  OC  n,  OC  J2,  and  OC  J3  are  large  enough,  due  to  (5) 
being  positive-definite,  the  minimization  is  stably  performed.  Thus, 
with  high  SNRs  reasonably  smoothed  2D  current  field  data  Jx  (i,j, 
K)  and  Jy  ( i,j ,  K)  are  obtained. 

Utilization  of  the  obtained  smooth  current  field  data  Jx 
(i,j,  K)  and  Jy  ( i,j ,  K)  and  their  2D  spatial  derivatives  calculated  by 
applying  the  differential  filter  [23]  with  an  appropriate  cut-off 
frequency  fc  to  the  current  filed  data,  high  SNR  approximations  of 
the  spatially  inhomogeneous  coefficients  in  (4)  are  obtained. 
Furthermore,  by  approximating  the  2D  gradient  operator  using  the 
first-order  finite  difference,  over  each  2D  ROI  £2  K  the  first-order 
finite-difference  equations  with  the  approximations  of  the 
inhomogeneous  coefficients  are  derived  with  exception  of  the  some 
boundary  points  of  £2  K  [24,  25].  Subsequently,  by  substituting  the 
reference  values  In  0  ;  into  the  finite-difference  equations 
representing  In  (7  ( i,j ,  K)  in  £2  K  crossing  CO  ,  (/  =  1,...,  N),  the 
following  system  of  equations  are  derived  to  determine  the 
unknown  2D  global  solution  s,  i,e.,  a  raw  vector  comprised  of 
unknown  globally  relative  conductivities  In  0  (i,j,K)  in  £2  K: 

JD’s  =  j,  (5) 

where  the  matrix  J  is  comprised  of  the  reasonably  smoothed  current 
densities,  the  matrix  D’  of  coefficients  of  the  first-order 
approximations  of  2D  partial  derivatives  of  In  0  (i,  j,  K),  and  the 
raw  vector  j  of  reasonably  smoothed  partial  derivatives  of  the 
current  densities.  Here,  we  should  keep  in  mind  that,  D’s  is 
designed  to  express  the  2D  locally  relative  conductivity  distribution 
[24], 


However,  the  current  field  data  and  their  spatial 
derivatives  in  J  and  j  has  remaining  estimation  errors.  Moreover, 
when  only  a  single  current  field  measurement  is  taken, 
configurations  of  current  sources  and  reference  regions  might  be 
improper  ones  [25].  That  is,  in  some  2D  ROIs  reference  regions 
crossing  CO  ,  (/  =  1 ,  ...  ,N)  may  be  short  in  the  parallel  direction  to 
the  imitated  2D  current  flow.  Such  configurations  lead  to 
inaccessible  regions  in  2D  ROIs;  thus  resulting  in  JD’  being 
numerically  singular.  These  current  estimation  errors  and  improper 
configurations  ate  coped  with  by  utilizing  the  a  priori  knowledge 
about  the  smoothness  of  the  target  distribution  In  0  (x,  y,  z), 

As  the  employed  gradient  operator  D’  has  a  unique 
least-squares  inverse  (D’tD’)  'D’t  with  respect  to  the  locally 
relative  2D  conductivity  distribution  D’s,  by  utilizing  the 
regularization  method,  we  just  only  make  the  distribution  D’s 
unique  and  smooth  enough  to  be  integrable  over  the  2D  ROI  £2  K. 
That  is,  in  addition  to  utilizing  the  fact  that  the  original  2D 
distribution  VlnCT  (x,y,Z)  is  bounded,  since  it  is  assumed  that 
V  In  0  (x,y,Z)  is  spatially  smooth,  we  obtain  the  following 
functional  with  respect  to  the  2D  global  distribution  s: 

Rs  (s)  =  llj  -  JD’sll2  +  Cf  sl  IID’sll2  +  as2IID’T  D’sll2,  (6) 
where  OC  sl  and  OC  s2  are  the  regularization  parameters  of  the 
penalty  terms,  i.e.,  L2-norms  of  the  gradient  and  Laplacian  of  s.  Due 
to  the  functional  (6)  derived  in  each  2D  ROI  £2  K  being 
positive-definite  even  though  some  2D  ROIs  have  originally 
inaccessible  regions,  numerical  minimization  of  the  every 
functional  yields  the  3D  globally  relative  conductivity  distribution 
through  stable  implicit-integration  over  an  arbitrary  path  of  the  2D 
locally  relative  conductivities  [23,  24]. 

Minimization  of  functionals  Rj  (J)  [(5)]  and  Rs  (s)  [(6)] 
are  performed  by  the  conjugate  gradient  method  [23,  24], 
Occasionally  the  respective  minimization  will  be  performed  in  3D 
ROI  £2  by  utilizing  the  approximation  of  3D  gradient  operator. 
The  suitability  of  this  combined  current  estimation  method  and 
conductivity  estimation  method  are  evaluated  in  the  next  section. 

III.  Evaluations 

Measured  tangential  magnetic  field  data  Bx  (i,j,  k0)  and  By  ( i,j ,  kg) 
are  simulated  by  current  flow  in  the  (x,  y)-plane  with  z  =  Z 
containing  a  central  circular  region  (dia.  =  10.0)  whose  conductivity 
value  differs  from  that  of  the  surrounding  medium,  i.e.,  4.00  versus 
1.00.  Specifically,  the  current  flow  is  generated  by  applying  current 
with  the  intensity  of  1.00  in  the  x-direction  at  infinite  distance,  and 
the  resultant  fields  Bx  (x,  y,  Zo)  and  By  (x,  y,  Zg)  are  supposedly 
measured  at  the  distance  from  0.2  to  1.0  over  a  finite  ROI  {  £2  K: 
-20  <  i  <  20,-20  <  j  <  20}  using  Ax=  Ay  =1.0.  The 
magnetic  fields  Bx  (i,  j,  k0)  and  By  O',  j,  kg)  are  filled  with 


independent  white  noises  of  the  same  power,  and  are  used  to 
estimate  the  current  density  fields  Jx  (i,  j,  K)  and  Jy  (i,  j,  K).  The 
SNRs  of  the  simulated  magnetic  field  data  By  (i,j,  K)  are  50. 40,  30, 
20,  and  10  dB.  Subsequently,  the  globally  relative  conductivity 
distribution  In  <7  (i,  j,  K)  in  £2  K  taken  to  have  a  proper 
configuration  of  a  reference  region  with  respect  to  the  applied 
current,  i.e.,  {-20  <  i  <  20,  7  =  -20}. 

In  the  respective  regularization  processes  in  (5)  and  (6), 
we  only  use  the  preferable  a  priori  knowledge  about  J  and  D’s. 
That  is,  CC  ji  and  CC  J2  in  (5),  and  CC  sl  in  (6)  are  always  set  equal 
to  zeros  though  CC  J3  and  CC  s2  are  varied.  In  the  conjugate  gradient 
method,  the  initial  estimates  of  J  and  s  are  set  as  a  spatial 
distributions  of  zeros. 

To  evaluate  the  accuracy  of  estimated  J  ( i ,  j,  K)  and 
In  <7  ( i,j ,  K),  we  use  SNRs  defined  as 
SNR|P(/J,  K)]  =  1 0  logio 

D  */,;,*)  I2 

_ _  (7) 

£  I  P(i,  j,  K)  -  P(i,  j,  K)  I2 

where  P  ( i ,  j,  K)  and  P  (i,  j,  K)  represent  the  calculated  and 

original  distributions,  respectively. 

Effectiveness  of  windowing  on  current  field  estimation 
is  shown  in  Fig.  1,  i.e.,  SNRs  with  no  application  of  the 
regularization  method  of  the  estimated  current  fields  Jx  (/,  j)  and  Jy 
( i,j )  vs.  the  distance  of  the  sensor  array  from  the  2D  ROI,  where  the 
SNR  of  the  no  windowed  magnetic  field  data  By  (i,  f)  is  a 
parameter:  00  dB  (raw  data),  40  dB,  and  lOdB.  Obviously, 
windowing  is  effective  particularly  when  the  SNRs  of  magnetic 
data  is  high.  The  distance  increasing,  the  accuracy  of  the  estimation 
significantly  degrades.  Figures.  2a  and  2b  respectively  show  Bx  (/, 
j)  (40dB)  and  By  (/',  j)  (16.1  dB)  when  the  distance  being  1.00, 
while  Fig.  2c  shows  windowed  By  ( i ,  j).  Figures.  2d  and  2e  are 
respectively  current  fields  Jx  ( i,j )  (20.0  dB)  and  Jy  (i,  j)  (-  0.597 
dB)  estimated  from  the  windowed  magnetic  field  data,  while  Fig.  2f 
shows  Jx  (i,  j)  estimated  from  no  windowed  magnetic  data. 

Effectiveness  of  regularization  on  current  field 
estimation  is  shown  in  Fig.  3,  i.e.,  regularization  parameter  values 
vs.  SNRs  of  the  estimated  current  fields.  The  estimates  Jx  (j,  j )  and 
Jy  (j,  j)  obtained  with  the  best  regularization  parameter  value  and 
with  no  regularization  are  shown  in  Figs.  2e-h,  6c-f,  7c-f,  and  8c-f 
respectively  when  the  SNRs  of  By  (i,j)  are  40,  30,  20,  and  10  dB. 

Effectiveness  of  regularization  on  conductivity 
estimation  is  shown  in  Figs.  4  and  5,  i.e.,  regularization  parameter 
values  vs.  SNRs  of  the  estimated  conductivity  distributions,  and 
SNR  of  By  (/,  j )  vs.  mean  values  of  the  best  estimated  conductivity 


values  at  the  central  squared  region  (3  x  3)  of  the  inhomogeneous 
circular  region.  The  best  reconstructions  0  (i,  j)  and  the 
reconstructions  with  no  regularization  are  shown  in  Figs.  2ij,  6gh, 
7gh,  and  8gh  respectively  when  the  SNRs  of  By  (i,j)  are  40,  30,  20, 
and  10  dB. 

IV.  Conclusions 

In  this  simulation,  to  quantitatively  evaluate  the  inhomogeneous 
conductivity  value,  synthetic  aperture  technique  should  realize  the 
accurate  measurements  of  tangential  magnetic  fields  with  SNRs 
larger  than  30  dB  [(7)]  though  the  inhomogeneous  circular  region 
could  be  detected  even  when  the  SNRs  of  the  magnetic  field  data 
are  very  low. 


Figure.  1 :  SNRs  of  the  estimated  current  densities  vs.  positioned  height  of 
magnetic  sensor  array. 


Figure.  2:  Current/conductivity  estimation  when  the  distance  being  1.0  and 
the  SNR  of  By  (i,j)  being  40  dB:  (a)  Bx  (i,j),  (b)  By  (i,j),  (c)  windowed  By 
(i,j),  (d)  Jx  (i,j)  and  (e)  Jy  (i,j)  obtained  from  windowed  magnetic  field  data 
with  no  regularization,  (f)  Jx  (i,j)  obtained  from  no  windowed  magnetic  data, 
(g)  Jx  (i,  j)  and  (h)  Jy  (j,  j)  obtained  with  regularization  from  windowed 
magnetic  data,  (i)  C 7  (i,  j )  obtained  with  no  regularization  from  the  current 
field  data  (g)  and  (h),  (j)  (7  (i,j)  obtained  with  regularization. 


Figure.  3:  Regularization  parameter  values  vs.  SNRs  of  the  estimated  current 
density  fields. 


Figure.  4:  Regularization  parameter  values  vs.  SNRs  of  the  estimated 
conductivity  distributions. 
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Figure.  5:  SNR  of  By  (i,j)  vs.  mean  value  of  the  best  estimated  conductivity 
values  at  central  squared  region  (3  x  3)  of  the  inhomogeneous  circular  region. 


Figure.  6:  Current/conductivity  estimation  from  windowed  magnetic  field 
data  when  the  SNR  of  By  (ij)  being  30  dB  (distance:  1.0):  windowed  (a)  Bx 
(i,  j )  and  (b)  By  (i,  j ),  (c)  Jx  (/,  j)  and  (d)  Jy  (/,  j)  obtained  with  no 
regularization,  (e)  Jx  (/,  j)  and  (f)  Jy  (i,  J)  obtained  with  regularization,  (g) 
(7  (i,  j )  obtained  with  no  regularization  from  current  field  data  (e)  and  (f), 
(h)  O"  (i,j)  obtained  with  regularization. 


Figure.  7 :  Current/conductivity  estimation,  the  SNR  of  By  (i,j)  being  20  dB. 


Figure.  8:  Current/conductivity  estimation,  the  SNR  of  By  (i,j)  being  10  dB. 
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